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ABSTRACT 

5 Equulei is among the most well-studied nearby binary star systems. Results of its observation 
have been applied to a wide range of fundamental studies of binary systems and stellar astrophysics. 
It is widely used to calibrate and constrain theoretical models of the physics of stars. We report 27 
high precision differential astrometry measurements of 5 Equulei from the Palomar High-precision 
Astrometric Search for Exoplanet Systems (PHASES). The median size of the minor axes of the 
uncertainty ellipses for these measurements is 26 micro-arcseconds (/ias). These data are combined 
with previously published radial velocity data and other previously published differential astrometry 
measurements using other techniques to produce a combined model for the system orbit. The distance 
to the system is determined to within a twentieth of a parsec and the component masses are determined 
at the level of a percent. The constraints on masses and distance are limited by the precisions of the 
radial velocity data; we outline plans improve this deficiency and discuss the outlook for further study 
of this binary. 

Subject headings: stars:individual((5 Equulei) - binaries: close - binaries: visual - tech- 

■ niques:interferometric - astrometry - stars:distances 
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■ The study of 5 Equulei as a binary star has lasted nearly 200 years. In the early 1800s, William Herschel's (mistaken) 
' listing of it as a wide binary (with what Friedrich Struve later proved to be an unrelated background star) brought 
, it to the attention of many astronomers. While making follow-up observations to support his father's claim that 

the proposed pair were only an optical double, Otto Wilhelm von Struve in 1852 found that while the separation of 
_ the optical double continued to grow (to 32"), the point-spread-function of 6 Equulei itself appeared elongated. He 
, concluded that 5 Equulei itself is a much more compact binary. (It is perhaps interesting to note that the Struve 
I ' family's study of the system continued through 1955, when Otto Struve and K. L. Franklin included the system in a 
O spectroscopic study.) 

5 Equulei (7 Equulei, HR 8123, HIP 104858, HD 202275, ADS 14773) is among the most well-studied nearby binary 



in 
o 
o 



>< 



1. INTRODUCTION 



^ star systems. It is particularly useful to studies of binary systems and stellar properties as it is close {d = 18.4 pc), 
bright {V — 4.49, K = 3.27), and can be studied both visually (semi-major axis roughly a quarter an arcsecond) 
^ . and spectros c opica lly (spectral classes F7V-I-F7V, K « 12.5 kms""'^) in a reasonable amount of time [P k, 5.7 years); 
iMazeh et al.l 1)199^ found only 23 binaries within 22 pc with periods less than 3000 days and spectra that were nearly 
solar (spanning types F7-G9, classes IV-V, V, and VI). As such, it is regularly included in statistical surveys of binary 
^ i systems (see, for example, Heacox 1998; Hale 19i).4) and fundamental stellar properties such as the mass-luminosity 
■ - ' relationship, calibrating photometric parallax scales, tabulat ing the H-R diagram of the solar neighborhood, and 
constraining models of stellar atmospheres (see, for example, 'Laste nn et et al.l l2002t lEggenl[l998t iCastelli et al.lll99^ 
ISmallcv and Dworctskv 1995; Bochm 1989; Habets and Hcintzc 19811 lPoDDeilll98n[r These applications depend upon 
accurate knowledge of the components' physical properties and the system's parallax. 

Previously, the visual orbit models (and thus evaluation of the total system mass and orbital parallax) of 5 Equulei 
have been limited by diffe rential astrometry with relativ e precisions on order of a few percent. The recently developed 
method by two of us fsee lLane and Muterspaughl 1)20041 ^) for ground-based differential astrometry at the 10 /ias level 
for sub-arcsecond ("speckle") binaries has been used to study 5 Equulei over the 2003-2004 observing seasons. These 
measurements represent an improvement in precision of over two orders of magnitude over previous work on this 
system. We present these new data, an updated three-dimensional model of the system, and the physical properties 
of the component stars. The observations were taken as part of the Palomar High-precision Astrometric Search for 
Exoplanet Systems (PHASES). The primary goal of PHASES is to search for planets in sub-arcsecond binaries, with 
the data also being useful for studies of binaries themselves. 

PHASES data is collected at the Palomar Testbed Interferometer (PTI). PTI is located on Palomar Mountain near 
San Diego, CA (Colavitact al. 1999). It was developed by the Jet Propulsion Laboratory, California Institute of 
Technology for NASA, as a testbed for interferometric techniques applicable to the Keck Interferometer and other 
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Division of Phy sics, Mathematics and Astronomy, 105-24, California Institut e of Technology, Pasadena, CA 91125 
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missions such as the Space Interferometry Mission, SIM. It operates in the J (1.2/xm), H (1.6/ini), and K (2.2/im) 
bands, and combines starUght from two out of three available 40-cm apertures. The apertures form a triangle with 
one 110 and two 87 meter baselines. 

2. OBSERVATIONS AND DATA PROCESSING 
2.1. PHASES Observations 

S Equulei was observed with PTI on 27 nights in 2003-2004 using the observing mode described in 
ILane a nd Mut erspaughl lj2004|) . We briefly review this method for phase-referenced differential astrometry of sub- 
arcsecond binaries here. 

In an optical interferometer light is collected at two or more apertures and brought to a central location where the 
beams are combined and a fringe pattern produced. For a broadband source of central wavelength A the fringe pattern 
is limited in extent and appears only when the optical paths through the arms of the interferometer are equalized 
to within a coherence length (A = A^/AA). For a two-aperture interferometer, neglecting dispersion, the intensity 
measured at one of the combined beams is given by 

/(:^)=/ofl + l/^^^^^sin(2WA)) (1) 

where x is the differential amount of path between arms of the interferometer, V is the fringe contrast or "visibility," 
which can be related to the morphology of the source, and AA is the optical bandwidth of the interferometer assuming 
a flat optical bandpass (for PTI AA = 0.4/im). 

The location of the resulting interference fringes are related to the position of the target star and the observing 
geometry via 

d^^ + 5a(y,t^ +c (2) 

where d is the optical path-length one must introduce between the two arms of the interferometer to find fringes. This 
quantity is often called the "delay." 1^ is the baseline, the vector connecting the two apertures. ~^ is the unit vector 
in the source direction, and c is a constant additional scalar delay introduced by the instrument. The term Sa , tj 

is related to the differential amount of path introduced by the atmosphere over each telescope due to variations in 
refractive index. For a 100-m baseline interferometer an astrometric precision of 10 /Lias corresponds to knowing d to 5 
nm, a difficult but not impossible proposition for all terms except that related to the atmospheric delay. Atmospheric 
turbulence, which changes over distances of tens of centimeters, angles on order tens of arcseconds, and on subsecond 
timescales, forces one to use very short exposures (to maintain fringe contrast) and hence limits the sensitivity of the 
instrument. It also severely limits the astrometric accuracy of a simple interferometer, at least over large sky-angles. 

However, in narrow-angle astrometry we are concerned with a close pair of stars, and we are attempting a differential 
astrometric measurement, i.e. we are interested in knowing the angle between the two stars (As — S2 — si). The 
atmospheric turbulence is correlated over small angles. If the measurements of the two stars are simultaneous, or 
nearly so, the atmospheric terms cancel. Hence it is still possible to obtain high precision "narrow-angle" astrometry. 

To correct for time-dependent fluctuations in the atmospheric turbulence, observations consisted of operating PTI 
in a phase-referenced observing mode. After movable mirrors in the beam-combining lab apply delay compensation to 
the light collected from two 40 cm apertures, the light from each aperture is split using 30/70 beamsplitters. Seventy 
percent of the light is sent to the phase-tracking "primary" interferometric beam combiner which measures the time- 
dependent phase of one star's interferogram (fringes) caused by the atmospheric turbulence, and used in a feed-forward 
loop to control the optical delay lines. 

The other 30% of the light is diverted to the "secondary" interferometric beam combiner. In this system we have an 
additional delay line with a travel of only w 500 microns. We use this to introduce delay with a sawtooth waveform 
with frequency on order a Hertz. This allows us to sample the interferograms of all stars in the one arcsecond detector 
field whose projected separations are within the scan range. Laser metrology is used along all starlight paths between 
the 30/70 split and the point of interferometric combination to monitor internal vibrations differential to the phase- 
referencing and scanning beam combiners. For 6 Equulei, the typical scanning rate in 2003 was one scan per second 
and four intensity measurements per ten milliseconds; these values were doubled in 2004. The typical scan amplitude 
was 100 microns. An average of 1700 scans were collected each night the star was observed over a time span of 30 to 
90 minutes. 

2.2. PHASES Data Reduction 

The quoted formal uncertainties in the PHASES data are derived using the standard PHASES data reduction 
algorithm, which we review here. First, detector calibrations (gain, bias, and background) are applied to the intensity 
measurements. Next, we construct a grid in differential right ascension and declination over which to search (in 
ICRS 2000.0 coordinates). For each point in the search grid we calculate the expected differential delay based on the 
interferometer location, baseline geometry, and time of observation for each scan. These conversions were simplified 
usi ng the routines from the Naval Observatory Vector Astrometry Subroutines C Language Version 2.0 (NOVAS-C; 
see lKaplan et al.l (^8^). A model of a double- fringe packet is then calculated and compared to the observed scan to 
derive a value as a merit of goodness-of-fit; this is repeated for each scan, co-adding all of the values associated 
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with that point in the search grid. The final surface as a function of differential right ascension and declination 
is thus derived. The best-fit astrometric position is found at the minimum position, with uncertainties defined by 
the appropriate contour — which depends on the number of degrees of freedom in the problem and the value of the 
X^-minimum. 

One potential complication with fitting a fringe to the data is that there are many local minima spaced at multiples 
of the operating wavelength. If one were to fit a fringe model to each scan separately and average (or fit an astro- 
metric model to) the resulting delays, one would be severely limited by this fringe ambiguity (for a 110-m baseline 
interferometer operating at 2.2/xm, the resulting positional ambiguity is ~ 4.1 milliarcseconds) . However, by using the 
X^-surface approach, and co-adding the probabilities associated with all possible delays for each scan, the ambiguity 
disappears. This is due to two things, the first being that co-adding simply improves the signal-to-noise ratio. Second, 
since the observations usually last for an hour or even longer, the associated baseline change due to Earth rotation 
also has the effect of "smearing" out all but the true global minimum. The final x^-surface does have dips separated 
by ~ 4.1 milliarcseconds from the true location, but any data sets for which these show up at the 4(T level are rejected. 
The final astrometry measurement and related uncertainties are derived by fitting only the 4ct region of the surface. 

The PHASES data reduction algorithm naturally accounts for contributions from photon and read-noise. Unmon- 
itored phase noise shows up by increasing the minimum value of x^ surface. Comparison of this value with that 
expected from the number of degrees of freedom allows us to co-add the phase noise to the fit. 

This method has been rigorously tested on both synthetic and real data. Data sets are divided into equal sized 
subsets which are analyzed separately. A Kolmogorov-Smirnov test shows the formal uncertainties from the PHASES 
data reduction pipeline to be consistent with the scatter between subsets. After an astrometric solution has been 
determined, one can revisit the individual scans and determine best-fit delay separations on a scan-by-scan basis 
(the fringe ambiguity now being removed). The differential delay residuals show normal (Gaussian )distribution, and 
Allan variances of delay residuals agree with the performance levels of the formal uncertainties and show the data to 
be uncorrelated. We conclude that the PHASES data reduction pipeline produces measurement uncertainties that 
are consistent on intranight timescales. Additionally, because the stars in b Equulei are nearly identical (and in 
particular have equal temperatures to within the uncertainties of the best measurements), potential systematics such 
as differential dispersion are negligible. 

The differential astrometry measurements are listed in Table in the ICRS 2000.0 reference frame. A Keplerian 
fit to the PTI data using the formal uncertainties found the minimized value of reduced Xr = 14.46, implying either 
that the uncertainty estimates are too low by a factor of 3.8 or the (single) Keplerian model is not appropriate for 
this system. Several possible sources of excess astrometric scatter have been evaluated. The PTI baseline vectors are 
determined by observing pointlike stars across the sky and inverting eq. |5| The baseline solutions are stable at the 
level of a millimeter over year-long timescales, a factor of ten better than is required for 10 microarcsecond astrometry 
on subarcsecond binaries. To test the effect of the interferogram model template on the astrometry, the data were 
refit several times using several models for the instrumental bandpass; the differential astrometry solutions varied by 
less than a tenth of the formal uncertainties. The maximum effect of starspots is evaluated as approximately 8 /xas 
of scatter per 10 millimagnitudes of photon ietric variability, a level n ot observed in b Equulei (Hipparcos photometry 
shows a scatter of only 4 millimagnitudes l)van Leeuwen et aDligQTj) '). The delay lines at PTI are in air rather than 
vacuum, introducing longitudinal dispersion to the system and color-dependent variations to the points of zero optical 
delay. Because the components of b Equulei are equal temperature, this effect cancels in a differential measurement. 
Because dispersion can couple to the astrometric measurement in more complicated ways (for example, if the optical 
alignment varies with either sky position or internal delay and the point spread function has color dependencies) and 
some PHASES targets do not have equal temperature components, we are implementing a dispersion compensator for 
operation beginning in the 2005 observing season. A more complete analysis of PHASES measurement uncertainties 
will be addressed in a future paper. 

The uncertainty values listed in Tabled have been increased by a factor of 3.8 over the formal uncertainties; these 
increased values are used in fits presented in this paper, in order that this data set can be combined with others. At 
this time we do not find that more complicated models (such as adding additional unseen system components) produce 
better fits to the PTI data. The rescaled (raw) median minor- and major-axis uncertainties are 26 (6.8) and 465 (122) 
/ias. The rescaled (raw) mean minor- and major-axis uncertainties are 35 (9.2) and 1116 (294) /ias. 

2.3. Previous Differential Astrometry Measurements 

We have collected previously publis hed differential astrom etry measurements made with other methods. Most of 
these measurements were tabulated bv lHartkopf et alJ l)2004f) in the Fourth Catalog of Interferometric Measurements 
of Binary Stars, though several additional measurements (particularly those made by micrometer measures) had to 
be researched from the original sources. In two cases we found di screpancies be t ween the uncertainties listed in the 
Fourth Catalog and the original sources (the 1977.8811 point bv iMorgan et alJ l(T98Qil and that from 1983.9305 by 
iBonneau et alJ l)1984H ): in each case we used the uncertainties listed in the original work. Several data points listed 
without uncertainty estimates in the Fourth Catalog were found to have uncertainty estimates listed in the original 
works, in which case those values were used. 

Most of the previous differential astrometry measurements were published without any associated uncertainties. 
To allow these to be used in combined fits with other data sets, we determined average uncertainties as follows. 
We separated the measurements into subgroups by observational method and analyzed each set individually; the 
first group included eyepiece and micrometer observations, and the second contained interferometric observations. 
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Table lU 
PHASES data for S Equulei 



JD-2400000.5 



(5RA 
(mas) 



(5Dec 
(mas) 



(^as) 



^ maj 

(/tas) 



(dcg) 



0"RA 
(/^as) 



""Dec 
(/ias) 



52896, 
52897, 
52915, 
52917, 
52918, 
52920, 
52929, 
52930, 
52951, 
53172, 
53173, 
53181, 
53182, 
53186, 
53187, 
53197, 
53198, 
53199, 
53208, 
53214, 
53215, 
53221, 
53222, 
53229, 
53236, 
53249, 
53251, 



18089 
15549 
20224 
18203 
17180 
19763 
15827 
16162 
10343 
44519 
40171 
39890 
39947 
39977 
39616 
39430 
36980 
39127 
31795 
33956 
33920 
26266 
34182 
30103 
23726 
24027 
19930 



43.1331 

42.9731 

46.9142 

46.8911 

46.9129 

46.8789 

49.5698 

49.3723 

53.1325 

30.9090 

30.6927 

28.0318 

27.4793 

25.5104 

26.2104 

21.9940 

21.8513 

22.2861 

17.7859 

15.9122 

15.4419 

14.0158 

12.6377 

10.0096 

8.2328 

1.7950 

2.4844 



28.3054 

29.2553 

41.1237 

42.6169 

43.3924 

44.7316 

50.3713 

51.1589 

64.3276 

105.4060 

104.8840 

101.4538 

101.1192 

99.5860 

98.5359 

93.8953 

93.2833 

92.4201 

87.9990 

84.5147 

83.9798 

80.9103 

79.9100 

75.6266 

70.6564 

62.9850 

60.8650 



27.9 

24.0 

34.8 

15.0 

26.1 

37.2 

18.0 

53.3 

119.7 

26.9 

34.4 

36.2 

29.3 

69.3 

22.7 

15.7 

14.9 

66.5 

19.7 

22.5 

13.6 

12.6 

16.3 

16.6 

97.6 

14.9 

54.3 



279.0 

378.4 

465.2 

273.3 

305.1 

1537.1 

333.9 

1344.3 

3391.8 

849.2 

629.2 

1586.5 

1339.1 

2784.6 

816.0 

126.8 

136.6 

2970.9 

197.4 

210.8 

132.7 

240.5 

784.0 

407.2 

3640.5 

192.5 

4775.1 



150.13 

146.47 

163.69 

160.60 

158.50 

166.02 

162.42 

163.85 

163.63 

152.99 

17.94 

149.30 

149.87 

151.54 

151.72 

155.88 

152.64 

156.61 

20.76 

154.02 

153.90 

15.60 

159.11 

154.97 

147.94 

153.47 

148.21 



242.3 

315.7 

446.6 

257.8 

284.1 

1491.6 

318.4 

1291.3 

3254.4 

756.7 

598.7 

1364.2 

1158.3 

2448.4 

718.7 

115.9 

121.5 

2727.0 

184.7 

189.8 

119.3 

231.7 

732.5 

369.0 

3085.8 

172.4 

4058.8 



141.1 

209.9 

134.8 

91.8 

114.4 

373.0 

102.3 

377.5 

962.9 

386.4 

196.6 

810.6 

672.6 

1328.2 

387.1 

53.8 

64.2 

1180.8 

72.3 

94.5 

59.6 

65.8 

279.9 

172.9 

1934.1 

87.0 

2516.1 



-0.97371 

-0.99054 

-0.96308 

-0.98499 

-0.96945 

-0.99472 

-0.98287 

-0.98912 

-0.99157 

-0.99695 

0.98292 

-0.99865 

-0.99873 

-0.99824 

-0.99778 

-0.94784 

-0.96511 

-0.99812 

0.95684 

-0.96443 

-0.96725 

0.98013 

-0.99806 

-0.99436 

-0.99822 

-0.98165 

-0.99968 



854 

1698 

885 

2197 

1181 

970 

1846 

584 

232 

1142 

1322 

679 

1391 

250 

2087 

2668 

3203 

515 

3244 

2792 

3947 

5719 

2302 

2316 

66 

2234 
794 



TABLE 1 

All quantities are in the ICRS 2000.0 reference frame. The uncertainty values presented in these data have all been 

SCALED BY A FACTOR OF 3.8 OVER THE FORMAL (INTERNAL) UNCERTAINTIES WITHIN EACH GIVEN NIGHT. COLUMN 6, <f>c, IS THE ANGLE 
between THE MAJOR AXIS OF THE UNCERTAINTY ELLIPSE AND THE RIGHT ASCENSION AXIS, MEASURED FROM INCREASING DIFFERENTIAL 
RIGHT ASCENSION THROUGH INCREASING DIFFERENTIAL DECLINATION (tHE POSITION ANGLE OF THE UNCERTAINTY ELLIPSE'S ORIENTATION IS 

90 — (jio)- The last column is the number of scans taken during a given night. The quadrant was chosen such that the 
larger fringe contrast is designated the primary (contrast is a combination of source luminosity and interferometric 

visibility). 



including speckle, phase-grating, aperture masking, and adaptive optics. The uncertainties were first estimated to be 
10 milliarcseconds in separation and 1 degree in position angle. A Keplerian model was fit to the data, and residuals 
in separation and position angle treated individually to update the estimates and outliers removed. This procedure 
was iterated until uncertainties were found consistent with the scatter. The 66 visual data points used have average 
uncertainties of 37.2 milliarcseconds in separation and 3.53 degrees. The 58 interferometric data points used have 
average uncertainties of 5.92 milliarcseconds and 1.59 degrees. 

We fit a Keplerian model to the data points for which uncertainty estimates were available to determine whether 
these were systematically too large or too small, and to find outliers. Because there were only four visual/micrometer 
measurements with published uncertainties, we did not treat these as a separate group. There were 42 interferometric 
measurements with published uncertainty estimates. We found the uncertainty estimates to be systematically too 
small; this factor was larger in position angle than in separation. Upon iteration, it was found that the separation 
uncertainties for these 46 data points needed to be increase by a factor of 1.71 and the position angle uncertainties by 
2.38. 

We list these previously published measurements in Tables [51 13 and 01 (complete versions available in the online 
version of this paper). 

2.4. Radial Velocity Data 

We have collected radial velocity data that h as previousl y been published by iDworetskv et al.l l)197lD and 
iPoDper and Dworctskv ( 1978) at Lick Observatory, Hans et alJ (|1979D at the Dominion Astrophysical Observatory 
(DAO), and Duquc nnov an d Mayor (1988) from CORAVEL. We list these measurements in Tables El El and |7| (full 
versions available in the online version of this paper). The Lick and DAO datasets were published without absolute 
uncertainty estimates but with relative weights assigned. We fit each data set independently to a Keplerian model and 
use the scatter in the residuals to determine the absolute uncertainties, which we present in the tables. A fit of the 
CORAVEL data to a Keplerian model showed excess scatter beyond the level of the formal uncertainties; we apply a 
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Table El 

Previous Astrometry With Uncertainties 



Year p 6 ap ag Reference 

(mas) (deg) (mas) (deg) 



1934.578 


161 


191.9 


8.55 


4.52 


1934.826 


100 


170.2 


17.10 


18.09 


1955.55 


380 


203.4 


17.10 


3.09 


1975.545 


117 


29.7 


30.78 


4.76 


1976.4496 


44 


247 


8.55 


3.57 


1976.6163 


90 


226.5 


7.69 


4.76 


1976.6217 


89 


225.9 


7.69 


4.76 


1977.6347 


293 


208.6 


7.52 


1.19 


1977.8811 


309 


202.8 


10.60 


4.76 




^Vilgm 
^ilson 
Worlev (19571 
Morgan et al. (1978'! 
McAlister (19781 
McAlister and Hendrv ('I982a1 
^cAlister^n^jiendr^ (1982a) 
McAlister and Hendrv d982a1 
Morgan et al. (19801 



TABLE 2 

Previous differential astrometry data with published uncertainties for <5 Equulei. p uncertainties have been increased by 
A factor of 1.71 and those for by a factor of 2.38. In many cas es 9 has been changed b y 180 degrees from the value 
appearing in the original works. Complete version available at [http: / / stuff.mit.edu /^^ atthewI / deltaEquTables / . 



Table El 

Previous Visual Astrometry Without Uncertainties 
Year p 8 Reference 
(mas) (deg) 



1852.63 370 196.3 Struve (18591 

1852.66 340 192.6 Struve (1859) 

1858.59 300 186.6 Struve (_1859) 

1858.59 300 188.9 Struve (18591 

1934.76 103 172 van den B os (19351 

1937.46 260 212.6 Wilson (1941V) 

1938.88 296 198.7 Danion (19521 

1939.513 262 200 Jeffers (b1 

1939.61 250 200.2 Jeffers (a) 

1939.690 250 199 Jcffcrs (b) 

TABLE 3 

Previous visual differential astrometry data without published uncertainties for S Equulei. The uncertainties presented 

IN THIS TABLE WERE DETERMINED BY THE SCATTER IN THE DATA. In MANY CASES 9 HAS BEEN CHANGED BY 180 DEGREES FROM THE VALUE 
APPEARING IN THE ORIGINAL WORKS. FOR THE VALUES IN THIS TABLE. ALL UNCERTAINTIES WERE TAKEN TO BE (Tp = 37.2 MAS AND 

= 3.53 DEGREES. Complete version available at |http://stuff.mit.edu/^^ atthew1/deltaEquTables/. 



Table El 

Previous Interferometric Astrometry Without Uncertainties 



Year 



Reference 



(mas) (deg) 



1979.5296 261 

1979.7700 227 

1980.4772 78 

1980.4798 74 

1981.4628 131 

1981.4710 132 

1981.7032 97 

1982.5059 136 

1982.7600 189 

1983.4232 292 



197.9 

195.4 

176.1 

175.3 

22.6 

21.2 

14.1 

217.3 

211.6 

206.7 



^c^iste^an^^iendr^ 

McAlister and Hendrv (1982b 1 



McAlister et al. 1 


'1983) 


McAlister ot al. 1 


1983) 


McAlister ot al. 1 


'19841 


McAlister et al. 1 


19841 


McAlist^e^^^ 1 


'19841 


McAlisto^^^ 1 


'1987a1 


McAlister et al. 1 


'1987a1 


McAlister et al. 1 


'1987a1 



TABLE 4 

Previous interferometric differential astrometry data without published uncertainties for 5 Equulei. The uncertainties 

PRESENTED IN THIS TABLE WERE DETERMINED BY THE SCATTER IN THE DATA. In MANY CASES 9 HAS BEEN CHANGED BY 180 DEGREES FROM 



AND ag = 1.59 DEGREES. COMPLETE VERSION AVAILABLE AT (hTTP: //STUFF. MIT. EDu/^^ ATTHEWl /dELTAEQuTABLES/ . 



scale factor of 1.527 to those uncertainties to allow these data to be combined with the other sets for simultaneous fits. 

3. ORBITAL MODELS 

The first correct orbital solution for 6 Equulei was that of [Luyte^ l)1934|1 . and is consistent with the modern 
orbit. Ivan de Kamp and LippincottI l)1945j) measured the astrometry of the binary photocenter and derived its first 
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Table m 

Lick Observatory Radial Velocities 



JD-2400000.5 RV 1 RV 2 

(kms~^) (kms^^) 



40410.5 


-6.8 


-24.7 


40428.3 


-5.2 


-24.9 


40430.3 


-4.9 


-26.0 


40459.3 


-3.4 


-27.3 


40487.2 


-1.6 


-28.8 


40494.3 


-1.7 


-29.4 


40516.2 


-0.8 


-31.4 


40544.1 


1.3 


-31.9 


40691.5 


-1.7 


-29.7 


40731.5 


-4.5 


-26.9 



TABLE 5 

Lick Observatory data from IDworetsky et al.I 1197111 and IPopper and D woretskyI 119781) used for the combined orbital fit 

IN THIS paper, a model FIT OF THIS DATA SET SEPARATELY WAS USED TO DETERMINE THE UNCERTAINTIES USED. NOTE THAT THE 
ORIGINAL WORK SWITCHED THE DESIGNATIONS OF THE PRIMARY AND SECONDARY COMPONENTS; THIS HAS BEEN CORRECTED HERE FOR 
COMBINATION WITH OTHER DATA SETS. FOR THE VALUES IN THIS TABLE, ALL UNCERTAINTIES WERE TAKEN TO BE 0.35 kms""^, EXCEPT FOR 

THE POINT ON JD — 2400000.5 = 4 0428.3 for which the unc ertainties were 0.70 kms"-*^. Complete version available at 

[HTTP : / / STUFF .MIT. EDU/ ~ [ MATTHEW 1 / DELTAEQuTABLES / . 



Table El 



DAO Radial Velocities 


Date 


RV 1 


o"i 


RV 2 




(years) 


(km s~ 


1) (kms- 


(kms~ 


1) (kms-1) 


1966.770 


-21.96 


0.81 


-9.70 


0.81 


1966.902 


-22.05 


0.81 


-9.73 


0.81 


1967.461 


-21.56 


0.81 


-10.17 


0.81 


1967.593 


-20.99 


0.81 


-9.87 


0.81 


1969.415 


-9.62 


0.81 


-22.18 


0.81 


1969.437 


-8.73 


0.81 


-23.01 


0.81 


1969.456 


-8.24 


0.41 


-23.69 


0.41 


1969.502 


-6.96 


0.41 


-24.48 


0.41 


1969.565 


-5.86 


0.41 


-26.12 


0.41 


1969.568 


-5.38 


0.41 


-25.46 


0.41 



T ABLE 6 

Dominion Astrophysical Observatory' data from IHans e t al. i 19791) used for the combined orbital fit in this paper. A 
model fit of this data set separately was used to determine the uncertainties used. Complete version available at 

HTTP: //stuff. MIT. EDU/~M ATTHEW1/dELTAEQuTABLES/. 



Table □ 
CORAVEL Radial Velocities 



JD-2400000.5 RV 1 cri RV 2 a2 

(kms^-"-) (kms^-^) (kms^-"-) (kms^-"-) 



43363.016 


-22.58 


0, 


.86 


-11.46 


0, 


.93 


43365.945 


-21.49 


0, 


.69 


-10.39 


0, 


.75 


43377.969 


-21.85 


0, 


.89 


-10.77 


1, 


.08 


43391.906 


-22.19 


0, 


.73 


-10.55 


0, 


.70 


43450.750 


-23.34 


1, 


.05 


-10.76 


1, 


.07 


43703.047 


-22.89 


0, 


.73 


-10.38 


0, 


.78 


43794.805 


-22.42 


1, 


.63 


-8.40 


1, 


.56 


43796.875 


-21.10 


0, 


.75 


-9.33 


0, 


.79 


43807.789 


-21.12 


0, 


.72 


-8.98 


0, 


.75 


44043.062 


-20.52 


0, 


.79 


-11.41 


0, 


.87 



TABLE 7 

CORAVEL data from IHuouennoy and IVIayorI 119881) used for the combined orbital fit in this PAPER. The uncertainties have 

BEEN REWEIGHTED BY A FACTOR OF 1.527 FROM THE ORIGINAL WORK, IN ORDER THAT THEY MIGHT BE COMBINED WITH OTHER DATA SETS 
FOR A SIMULTANEOUS FIT. NOTE THAT THE ORIGINAL WORK SWITCHED THE DESIGNATIONS OF THE PRIMARY AND SECONDARY COMPONENTS; 
THIS HAS BEEN CORRECTED HERE FOR COMBINATION WITH OTHER DATA SETS. COMPLETE VERSION AVAILABLE AT 
|HTTP: / /stuff. MIT. EDU/ ~^ ATTHEWl /dELTAEQuTABLES/. 
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Tablcigi 

Orbit models for S Equulci. 





C: -'-A ^ j-1 T It ■: n 1 >-n / 1 Q Q Q \ 

CI c? r 11 1 c^^n ^^^^^yj 










irriAjiijj -f- 












-j- RV 


Pre. ± RV 


-P (yr) 


5.713 


5.703 ±0.0070 


5.706 ±0.003 


5.40 ±0.11 


5.7059 ±0.0003 


5.7058 ±0.0003 


P (days) 






2084.08 ±0.92 


1974 ±39 


2084.07 ±0.12 


2084.05 ±0.11 


To (yr) 


1992.85 


1981.47 ±0.012 


2004.285 ±0.015 


2004.299 ±0.002 


2004.2954 ±0.001 


2004.2950 ±0.001 


To (MJD) 






53109.9 ±5.5 


53114.53 ±0.75 


53112.90 ±0.45 


53112.76 ±0.42 




0.44 


0.440 ±0.0046 


0.4519 ±0.0029 


0.415 ±0.0080 


0.437001 ±0.000076 


0.436983 ±0.000072 


a (mas) 


231 


232 ±1.8 




222.8 ±3.2 


231.85 ±0.11 


231.88 ±0.11 


^0,Liak (kms-1) 






-15.398 ±0.097 




-15.40 ±0.10 


-15.40 ±0.11 


vi 1 L/ j4. L-^ ■ ' 






-15.876 ±0.074 




-15.875 ±0.078 


-15.875 ±0.081 


^ (km s~ 1 




-15.85 ±0.074 


-15.728 ±0.095 




-15.73 ±0.10 


-15. 73 ±0.10 


Ml ( Mr~. 1 




1.19 ±0.034 






1.192 ±0.012 


1.193 ±0.012 


JW 2 W^f J 




1 12 ±0 032 






1 1 s"? 4-n ni 1 

-1- . i- O / JIU . KJ i- i- 


1 188 ±0 012 


Ml + M2 (Mq) 


2.35 








2.380 ±0.018 


2.380 ±0.019 


M1/M2 


1.06 ±0.03 


1.06 ±0.018 


1.004 ±0.011 




1.004 ±0.012 


1.004 ±0.012 


i (dcg) 


99 


99.0 ±0.43 




99.520 ±0.052 


99.394 ±0.020 


99.396 ±0.019 


^ (dog) 


3 


8.0 ±1.0 


187.01 ±0.91 


188.53 ±0.25 


187.96 ±0.12 


187.92 ±0.11 


f2 (deg) 


203 


203.8 ±0.29 




203.292 ±0.044 


203.301 ±0.046 


203.312 ±0.043 


TT (mas) 


54.32 ±0.90 


55.0 ±0.67 






54.38 ±0.14 


54.39 ±0.15 


d (pc) 










18.388 ±0.048 


18.386 ±0.050 


a (AU) 










4.263 (Derived) 


4.263 (Derived) 


asini (AU) 






4.186 ±0.011 




4.206 (Derived) 


4.206 (Derived) 


i^l (kms-1) 






12.224 ±0.074 




12.181 (Derived) 


12.182 (Derived) 


iC2 (kms-1) 






12.272 ±0.074 




12.229 (Derived) 


12.230 (Derived) 



TABLE 8 

Pre.: Previous differential astrometry measurements, listed in Tables|^|21 and0 

In the combined PHASES + previous + RV fit, all parameter uncertainties have been increased by a factor of y^Xr = \/l-17 

(though the Xr of the combined fit is artificial due to rescaling the uncertainties of the individual data sets, this 
reflects the degree to which the data sets agree with eachother). 



photocentric orbit. Their measurements yielded a measure of the mass ratio of 0.508:0.492 ±0.016 l|van de Kamd 
((l954i|) later also derived individual masses of 1.96 and 1.89 Mq with the same method, values which are too large due 
to an underestimated parallax). The first spectroscopic orbit was bv lDworetskv et al.l (119711). providing a ma ss ratio 
of roughly 1.044. Finally, a full three-dimensional model for the system was determined bvPHans et alTlll979l). Since 
that time, several more orbital solutions have been offered (see, for example. lStarikov'alll981t iDuauennov and Mavoil 
IT988llHa,rtkonf et a,].l [T996l lSoderhie]mllT999tlPourba ixl l200?l . 

We introduce three model parameters for the system velocity Vo, one corresponding to each obse rvatory frorn which 
radial velocity data is obtained. We do this to allow for instrumental variations; in particular, 'H ans et alJ l)1979f) 
notes a potential zero-point discrepancy of 500 ms~^ in data sets. Having fit each data set independently to correct 
uncertainty estimates, we combined all into a simultaneous fit to best determine system parameters. The results are 
listed in Table |H1 and plotted in Figure Q 

Each fit was repeated several times varying the set of non-degenerate parameters used in order to obtain uncertainty 
estimates for a number of desired quantities. The fit to radial velocity data alone was fit once using asini and 
R = M1/M2 as parameters, and again replacing these with Ki and K2, the velocity amplitudes. Similarly, the 
combined fits were repeated replacing parameters {a, R} with the sets {M = Mi + M2, R} and {Mi, M2}. Quantities 
in the combined fit that were derived from other parameters are listed separately at the end of Table |H1 

Despite spanning less than a year (a sixth of the orbit), the PHASES data by themselves are able to constrain 
many orbital parameters better than previous observations. We note in particular that the orbital angles are very well 
constrained. However, the relatively short time coverage of the PHASES data presents strong degeneracies between 
the system period, eccentricity, and semi-major axis, which increases the fit uncertainties to levels much larger than 
one would expect given the precision of the astrometry. If, for example, one holds fixed the period and eccentricity at 
the fit values, the uncertainty in semi-major axis drops from 3200 to 92 /ias. 

It is noted that the PHASES only fit results in values of period, eccentricity, and semimajor axis that agree with 
previous fits at only the 3ct level. At this time the causes of these discrepancies are not yet known. As noted, these 
parameters are degenerate with each other in the PHASES only fit, thus it is not surprising that all are discrepant 
at the same level. However, the data sets agree well enough to be simultaneously fit in a combined model (for which 
parameter degeneracies are less of a problem), and the resulting parameters do agree well with previous orbital models. 

We find that the addition of the previous differential astrometry to the combined model does little to improve the 
fit. The combined fit is limited by the precision of the radial velocity observations. Using the techni que for ob taining 
high precision radial velocimetry on double line spectroscopic binaries using an iodine cell reported in lKonackil (|2P04) , 
one of us (MK) is leading an effort to obtain such data to better constrain the orbital model, component masses, 
and system distance. The PHASES program will also continue to observe S Equulei so that the combination of high 
precision radial velocities and differential astrometry can be used for a comprehensive search for giant planets orbiting 
either star. We have conducted simulations which show collecting ten radial velocity measurements with 20 ms~^ 
precisions during the 2005 observing season will improve the constraints on the component masses by a factor of two. 
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Fig. 1. — The orbit of <5 Equulei. In the upper left plot, only interferometric differential astrometry measurements for which the original 
works offered measurement uncertainty estimates have been included. Note also that the dimensions of these uncertainty ellipses are those 
quoted by the original authors. The dimensions of the uncertainty ellipses plotted for the PHASES measurements have been stretched 
by a factor of 3.8 as discussed in the text. The high ellipticities of PHASES uncertainty ellipses are caused by use of a single baseline 
interferometer coupled with the limited range of hour angles over which 5 Equulei could be observed (due to limited optical delay range at 
PTI). 



4. PARALLAX 

Early attempts to measure the parallax of b Equulei suffered from sy stematic errors due to its bina ry nature until 
iLuytenI determined a model for the visual orbit. In the same paper that lyan de Kamp and Lippincot^ determined the 
first orbit of the (5 Equulei photocenter, they determined a trigonometric parallax of 48 ifc 5 mill i arcsec onds. 

The best current values for the trigonometric parallax of (5 Equulei are given by IGatewool l|1994l) as an average 
of ground based observations (54.2 ± 0.93 mill-arcseconds) and from the Hipparcos mission (a binary-orbit corrected 
parallax of 54.32 ± 0.90 milliarcseconds is reported (Soderhielm 1999)). Our combined orbital solution using PHASES 
differential astrometry and previously published radial velocity measurements provides the best estimate of orbital 
parallax, at 54.39 ± 0.15 milliarcseconds, in good agreemen t with the trigonometric values and the previous best 
orbital parallax of 55 ± 0.67 milliarcseconds l|Pourbai3dl2000() . The orbital parallax determination is limited by the 
precision of the radial velocity measurements; simulations show that the high precision radial velocity observations 
planned for the next observing season will improve the precision by a factor of two. 

5. SYSTEM AGE 

The measured apparent V magnitude for the system is reported as 4.487±0.02 by the S imbad astronomical database . 
Combining this value with the measured ISV between the stars of 0.09±0.04 measured bv lten Brummelaar et alJ l)2000f) 
using adaptive optics on the Mt. Wilson 100" telescope and the distance determined by our orbital model, we find the 
components have absolute magnitudes o f V^ = 3.87 ±0.028 and V2 — 3.96 ±0.029. These values are combined with the 
stellar evolution m odels o flGirard i et al.l |2000) to determi ne the system age. The system's metallicities (in solar units ) 
of [Z/H] = -0.07 ()Grav et al.ll2003() and [Fe/H] = -0.07 l|Nordstrom et alJl2004j) most closely match iGirardi et alf s 
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Fig. 2. — Isochrones for stars of near-solar metallicities as functions of stellar mass and absolute magnitude. Isochrones shown are 
separated by logt = 0.05. Also plotted are the properties of the two components of 6 Equulei. The system age is logi = 9.35^q :[g. 



isochrone for stars of solar metallicity. We find the system age is logi = 9.35lg;i5 (~ 2.2 ± 0.6 Gyr). The relevant 
isochrones are plotted in Figure El 



S Equul ei is a sample system disc usse d in the S-type (orbiting just one stellar component of a binar y) planet stability 
studies of iRabl and DvorakI lIlQS S') and'Holm an and Wiegern l|199fllh . The numerical simulations of iR.abl and Dvorak! 
determined that planets were stable around either star if their orbital semi-major axis were 0. 68 AU (P = .34 year) 
or smaller; an additional semi-stable region existed out to 0.86 AU. The conclusion of Holmai Tand Wiegerti was that 
the regions of stability were of size 0.67 AU {P = 0.43 years; they assumed slightly different values for the component 
masses) around the primary and 0.66 AU around the secondary {P ~ 0.42 years). From these studies, we assume a 
stable region of roughly | AU around each star in which planets could be found. 

We do not find any obvious periodicities in our residuals, which are plotted in Figures 13 El and|Sl Periodograms of 
the PHASES fit residuals show no clear peaks between one and 180 days. We attempted to refit the PHASES data 
using a double-Keplerian model; we started each attempt with seed values for the wide Keplerian portion equal to the 
values found for the single Keplerian fit, and chose 3500 different starting values between one and 200 days for the 
period of the narrow portion (secondary Keplerian). The final value for the reduced Xr never found to fall below 
12.69, which is not significantly different from the value from the single Keplerian model of 14.46; because several seed 
periodicities produced Xr near the 12.69 level, we conclude this slight dip is a result of random noise and the data 
sampling function. We can safely conclude that there are no periodic signals in our residuals at the level of 100/xas, 
at least along the average minor axis of our uncertainty ellipses. At this conservative level we can conclude that there 
are not additional companions of mass 



The orbit of a third body could be hidden if it happens to be high inclination and coaligned with the major axis of our 
uncertainty ellipses. A more thorough analysis of the fit residuals and better constraints on companion masses will be 
part of a future investigation. 



The high precision differential astrometry measurements of the PHASES program are used to constrain the distance 
to 6 Equulei more than four times more precisely than previous studies, despite covering only a sixth of the orbit. The 
orbital parallax agrees well with the trigonometric one determined by Hipparcos observations. Whereas characterization 
of the system was previously limited by the precisions of differential astrometry measurements, it is now limited by 
the radial velocity observations. Continued monitoring of this nearby standard binary will be useful to search for 
additional system components as small as a Jupiter mass in dynamically stable orbits. 

We thank the support of the PTI collaboration, whose members have contributed designing an extremely reliable 
instrument for obtaining precision astrometric measurements. We acknowledge the extraordinary efforts of K. Rykoski, 
whose work in operating and maintaining PTI is invaluable and goes far beyond the call of duty. Observations with PTI 
are made possible through the efforts of the PTI Collaboration, which we acknowledge. Part of the work described in 
this paper was performed at the Jet Propulsion Laboratory under contract with the National Aeronautics and Space 
Administration. Interferometer data was obtained at the Palomar Observatory using the NASA Palomar Testbed 
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Fig. 3. — Residuals for PHASES differential astrometry of S Equulei. The error bars plotted have been stretched by a factor of 3.8 over 
the formal uncertainties as discussed in the text. The high cUipticity of the uncertainty ellipses causes neither the right ascension nor the 
declination uncertainties to be near the precision of the minor axis uncertainties, which have median uncertainty of 26 ^as. Due to the 
roughly North-South alignment of the baseline used for 24 of the 27 measurements, our more sensitive axis was typically declination. The 
bottom left plot shows the residuals along a direction that is 154 degrees from increasing differential right ascension through increasing 
differential declination (equivalent to position angle 296 degrees), which corresponds to the median direction of the minor axis of the 
PHASES uncertainty ellipses. Because the orientation of the PHASES uncertainty ellipses varies from night to night, no single axis is ideal 
for exhibiting the PHASES precisions, but this median axis is best aligned to do so. The bottom right plot shows residuals along the minor 
axis of each measurement's uncertainty ellipse. 
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Fig. 4. — Residuals for previous differential astrometry of <5 Equulei. 
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